Symmetry of carrier-envelope phase difference effects in strong-field, few-cycle 

ionization of atoms and molecules 
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In few-cycle pulses, the exact value of the carrier-envelope phase difference (CEPD) has a pro- 
nounced influence on the ionization dynamics of atoms and molecules. We show that for atoms in 
circularly polarized light, a change in the CEPD is mapped uniquely to an overall rotation of the 
system, and results for arbitrary CEPD are obtained by rotation of the results from a single calcu- 
lation with fixed CEPD. For molecules this is true only for linear molecules aligned parallel with the 
propagation direction of the field. The effects of CEPD are classified as geometric or non-geometric. 
The observations are exemplified by strong-field calculations on hydrogen. 
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Nowadays, it is possible to construct and control in- 
tense laser pulses with only a few optical cycles [l|, 
i.e., pulses described by a vector potential of the form 
A(t) = A f(t)sin(u)(t — J) + <ft)e, where A is the am- 
plitude, f(t) is the envelope, r is the pulse length, u is 
the frequency, <f> is the carrier-envelope phase difference 
(CEPD), and e is the polarization vector. The corre- 
sponding electric field is obtained from the vector poten- 
tial by E{t) = —dtA(t), and is shown in Fig. [T]for a sin 2 
envelope. Such pulses can be used to probe molecular and 
atomic dynamics on a very short time scale @, Q. The 
associated ionization dynamics becomes sensitive to the 
exact shape of the pulse and the carrier-envelope phase 
difference (CEPD) (see Fig. [T]). This dependence may 
be understood by the exponential dependence of the ion- 
ization rate on the instantaneous field strength and the 
corresponding emergence of the electron into the field- 
dressed continuum at specific instants of time during the 
pulse 0]. Asymmetries in the photoelectron spectrum 
may give information about the CEPD, and hence help 
in the characterization of the field. Electrons released at 
different times are accelerated to different final momenta 
and this fact is exploited in attosecond streaking [f| to 
map the time distribution of the pulse into a momen- 
tum distribution of the photoelectron, and to charac- 
terize the ultra-fast pulses [fj. Hence, just as few-cycle 
pulses are diagnostic tools for atoms and small molecules, 
the very same systems serve as diagnostic tools for the 
pulses themselves 0- The latter statement, of course, as- 
sumes that an accurate theoretical description is at hand 
for the pulsed laser-matter interaction. It is the purpose 
of this work to add further to this understanding. In par- 
ticular we are concerned with the CEPD effects, and a 
geometric interpretation of these. 

Carrier-envelope phase difference effects were stud- 
ied theoretically with emphasis on CEPD-induced spa- 
tial asymmetries in the ATI-spectrum/angular distribu- 
tion in a number of pap ers on strong- field ionization of 
atoms [1, 0, fiil U , Il2l ] and molecules [3] , strong-field 
dissociation |l4l ] , and high- harmonic generation fla j . The 



FIG. f : (Color online). The electric field E(t) /Eo, normalized 
to the peak field strength Eo and the pulse envelope f(t) as a 
function of time for two values of the CEPD, <f>. The electric 
field points in opposite directions for = and <f> = it. In ion- 
ization, the electric field will shake the electron until it gains 
enough energy to escape the Coulomb potential, and the an- 
gular distribution will depend on CEPD because the electric 
field (and the force F = —E) points in opposite directions for 
cj> — and <j> = 7T. The field parameters are r = 441 a.u. and 
lu = 0.057 a.u. (800 nm). 



asymmetries can be used to extract information about 
and ultimately to measure the CEPD. In Ref. jl6[ a spa- 
tial asymmetry in ionization with few-cycle circular po- 
larized laser pidses was observed for the first time. In 
Ref. 



17j | the generation of intense few-cycle laser pulses 



with stable CEPD was demonstrated, and a way of mea- 
suring CEPD, based on soft-X-ray radiation was pre- 
sented. In Ref. [l8| a spatial asymmetry in the ionization 
with few-cycle linear polarized laser pulses was measured, 
and the CEPD was determined with an estimated error 
of 7r/10. Asymmetries in the ionization signal combined 
with an attosecond pump pulse @ were used to measure 
directly the field of a linearly polarized few-cycle pulse 0] • 
Very recently, the detailed control over the CEPD played 
a crucial role in attosecond electron dynamics (IR Hfjl. lifj] . 
and experiments supporting the present findings have 
been reported at conferences [ilj ]. 

In this work, we present a systematic theoretical study 
of the symmetry of the response of atomic and molecular 
systems under a change of the CEPD in few-cycle pulses, 
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and we exemplify the discussion with a study of the dif- 
ferential electron momentum distribution for hydrogen 
under such pulses. 

We start out by considering an n-electron atom inter- 
acting with a few-cycle circularly polarized laser pulse 
described by the vector potential 



(1) 



x ^cos(?7 + 6 - 



^-)e x + sin(ry + . 



with f(rj) = sin (^y) the envelope, N the number of 
optical cycles, r\ = ujt — kz, uj the frequency, and k = ke z 
the wave vector. In the present case, with full inclusion of 
the spatial dependence of the field, the interaction of the 
atom with the field is obtained by the minimal coupling 
p — > p + A and the time-dependent Schrodinger equation 
reads [atomic units (a.u.) with m e — e = ao = h — 1 are 
used throughout], 
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where Hq is the field-free Hamiltonian. We are interested 
in relating this equation to an equation for the 6 = 
case. To this end, we note that the unitary operator 
D(6) = exp(jj z 0), where J z is the total angular mo- 
mentum, corresponds to a rotation of our system by an 
angle —6 around the z-axis. Since Hq is invariant un- 
der rotations around the z-axis, the transformed wave 
function = D{4>)'^ satisfies the Shrodinger equation, 

i-V =(ffo + J>(0)^,f,fi) -ffj^{4>) (3) 



nA 



^D(0)/ 2 (^)I>t(0))*', 



and using the Baker-Hausdorff lemma [22j we obtain, 
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^ = | H + ^ A(o, t, fj) ■ ft + ) 

(4) 

When we compare @ and Q , we see that a change in the 
CEPD from <f> — to 6 = 6 ' corresponds to a rotation 
of our system around the z-axis by the angle 6 ' ' ■ This 
fact is not only theoretically interesting, but also helpful 
in practical calculations since this symmetry property re- 
duces the number of computations one has to perform to 
a single one - all other results are obtained by suitable 
rotations. For instance, imagine we are interested in the 
differential ionization probability ^ ^ (<f>) for the mo- 
menta q x , q y in the polarization plane and for a general 
4> > 0. Then we calculate ^qdq = an< ^ r °t a te the 
result counterclockwise by an angle (f> to obtain d ^ P ^ q {4>) ■ 
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In the above derivation it is essential that the field-free 
Hamiltonian Hq is invariant to rotations around the z- 
axis. If this is not the case, the proof breaks down. As 
an example we look at ionization of diatomic molecules, 
or more generally linear molecules, in circularly polar- 
ized few-cycle laser pulses. If the molecule is aligned 
along the laser propagation direction, then the field-free 
Hamiltonian still has the required symmetry and the the- 
orem holds. If the molecule is not aligned along this 
axis, the system does not have the required symmetry 
and CEPD effects can not be reduced to a geometrical 
rotation. Accordingly, we may make a distinction be- 
tween rotational invariant atomic and molecular systems 
where the CEPD effects are purely geometric rotations, 
and systems which are not rotational invariant in which 
case true non- geometrical CEPD effects occur. For ex- 
ample, the results on strong-field ionization of with 
the molecule in the polarization plane [l3| belong to the 
latter category. 

As an illustration of the present findings, we consider 
ionization of atomic hydrogen, H(ls), in the strong-field 
approximation (SFA) [23|. We assume the dipole ap- 
proximation which means that rj — tot in (fTJ) and use 
that E — —d t A. The probability amplitude for direct 
ionization reads 



T fi = -i 



^ff{r,t)\E-r\Mr,t))dt, 



(5) 



where ^(r, t) is the Is ground state wave function. The 
final state f(r,t) is represented by a Coulomb- Volkov 
wave function [24j with asymptotic momentum q. The in- 
tegral in |5]| is analyzed as in [24[ and the numerical eval- 
uation is performed with Gauss-Legendre quadrature. 
By exploiting the invariance of the dot-product under ro- 
tations (E{4>) ■r*= [R z (<p)E(<p = 0)] • [R z (^)R z {-^)r\ = 
E(4> = 0) • f ', where R z (4>) is the 3x3 matrix describing 
rotations around the z axis), the SFA is readily shown to 
respond to CEPD changes like the exact theory discussed 
above. 

Once Tfi is known, a simple numerical integration over 
q z gives us the (q Xl q y ) distribution ^'^ = J\T fi \ 2 dq z . 
Figure presents the calculated distribution for various 
values of cf>. For varying cj>, the distribution rotates in 
accordance with the general theory. For 6 — 0, f Fi , % 

(IQ.r o Q u 

has a peak around (q x = 0, q y ~ 0.6) (this peak corre- 
sponds to about 14 photons above threshold), and is al- 
most symmetric around the line q x = 0. As 6 increases, 
this line is turned counterclockwise, so that has a 

peak around (q y = 0, q x ~ —0.6) for 6 = 

Some of the features can be explained by a semiclas- 
sical two-step model llj: first the electron escapes to 
the continuum at t = to with velocity v(to) = 0. Sec- 
ond, it moves like a classical particle under the influ- 
ence of the external Coulomb and laser fields. If, we 
neglect the Coulomb potential, which is justified if the 
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FIG. 4: (Color online) The {q x ,q y ) distribution for strong- 
field ionization of H(ls) by a linearly polarized field for (j> — 



0,7r/2, with I = 5.0 x 10 13 W/cm 2 



0.057 corresponding 



to a central wavelength of 800nm and three cycles, iV 



The grid size is dq x = dq y 



0.02. 



FIG. 2: (Color online) The (q x ,q y ) distribution for strong- 
field ionization of H(ls) for various values of <j>, with I = 5.0 x 
10 13 W/cm 2 , to — 0.057 corresponding to a central wavelength 
of 800nm and three cycles, N = 3. The grid size is dq x = 
dq v = 0.01. 
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FIG. 3: The curves show the time-dependent vector potential 
and the lines show if* for cj> = 0, tt/2. Laser parameters are as 
in Fig. 



field is very strong, Newton's second law tells us that 
q x (t) = A x (t)+q d , q y {t) = A y (t)+q d where q d x = -A x {t ) 
and q d = — A y (to). This means that the photoelectron 
direction, after the pulse where A vanishes, is determined 
by the vector q d = (q d , q d ). For simplicity t is chosen to 

be t/2, which is the time when \E(t)\ reaches its maxi- 
mum. Figure [3] shows plots of the time-dependent vector 
potential and q d = —A(to) for various values of <j>. The 
model explains how the preferred direction(the line of 
symmetry) of the photoelectron depends on CEPD. 

Finally we consider the case of a linearly polarized laser 
field and return to the question of symmetry. The vector 
potential in the linear case, A((j),t,r) — A f(r]) sin(r; + 
+ %/2)e y , can be written as A((f>, t, r) — t) — 

A R (^,t)), where 1 L/R (<f>,t,r) = ^f( v )(cos(i 1 + + 
ir /2)e x ±sm(r] + (j)+Tr /2)e y ) describe left(+) and right (—) 
hand circularly polarized fields. This relation suggests 
that we can characterize CEPD effects for a linearly po- 



larized pulse similarly to circularly polarized pulses. Un- 
fortunately, in exact theory, it is not possible to separate 
the effects of the circularly polarized fields, since they 
do not individually commute with the field-free Hamilto- 
nian. 

Figure 4 shows two examples for a linear field correspond- 
ing to the 4> = and <f> = tt/2 cases. These results may 
be interpreted in terms of the classical model discussed 
above (q y — ~A y {to)). For <j) — 0, the vector potential 
peaks in the negative y direction at the peak of the en- 
velope and produces electrons with positive momentum, 
q y = —A y (^). For <p = 7r /2, the vector potential is anti- 
symmetric with respect to the peak of the envelope, and 
the electron spectrum reflects this in a forward-backward 
symmetry. For <f> ~ 7r/2, the spectrum has moved to 
smaller momenta since the two maxima of the vector po- 
tential are lower in this case. Notice, however, that the 
model can not explain the size of the shift. There are 
two reasons for this discrepancy. First, the model ne- 
glects the Coulomb potential. Second we assume that 
vo = 0, which is not necessarily correct. 

In this work, we provided an exact characterization 
of CEPD effects in the interaction between circularly 
polarized few-cycle pulses and multielectron atoms and 
molecules including full account of the spatial depen- 
dence of the field (non-dipolc effects). For systems which 
are invariant to rotations around the propagation direc- 
tion of the field, a change of the CEPD corresponds to 
a rotation of the systems. For rotationally non-invariant 
systems, true non-geometrical CEPD effects may be ob- 
served. Calculations on H(ls) quantified the effects, and 
showed that observation of the momentum distribution 
can be used to measure CEPD for a circularly polarized 
laser pulse. For linearly polarized fields, the electrons 
spectrum was interpreted in terms of a classical model. 

There is currently considerable interest in control 
of the response of atoms and small molecules to a 
strong few-cycle pulse [H, e.g., in connection with align- 
ment dependent high-harmonic generation from small 
molecules [2|| [26[. Here the harmonics provide a deli- 
cate probe of the molecular wave packet dynamics and, 
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in turn, steering the electronic wave packets controls the 
harmonics. The carrier-envelope-phasc difference is one 
of the essential control parameters in few-cycle pulse and 
attosecond science. 

We thank T.K. Kjeldsen and U.V. Poulscn for useful 
discussion. L.B.M. was supported by the Danish Natural 
Science Research Council (Grant No. 21-03-0163) and 
the Danish Research Agency (Grant. No. 2117-05-0081). 
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Nowadays, it is possible to construct and control in- 
tense laser pulses with only a few optical cycles [l|, 
i.e., pulses described by a vector potential of the form 
A(t) = A f(t)sin(u)(t — J) + <ft)e, where A is the am- 
plitude, f(t) is the envelope, r is the pulse length, u is 
the frequency, <f> is the carrier-envelope phase difference 
(CEPD), and e is the polarization vector. The corre- 
sponding electric field is obtained from the vector poten- 
tial by E{t) = —dtA(t), and is shown in Fig. [T]for a sin 2 
envelope. Such pulses can be used to probe molecular and 
atomic dynamics on a very short time scale @, Q. The 
associated ionization dynamics becomes sensitive to the 
exact shape of the pulse and the carrier-envelope phase 
difference (CEPD) (see Fig. [T]). This dependence may 
be understood by the exponential dependence of the ion- 
ization rate on the instantaneous field strength and the 
corresponding emergence of the electron into the field- 
dressed continuum at specific instants of time during the 
pulse 0]. Asymmetries in the photoelectron spectrum 
may give information about the CEPD, and hence help 
in the characterization of the field. Electrons released at 
different times are accelerated to different final momenta 
and this fact is exploited in attosecond streaking [f| to 
map the time distribution of the pulse into a momen- 
tum distribution of the photoelectron, and to charac- 
terize the ultra-fast pulses [fj. Hence, just as few-cycle 
pulses are diagnostic tools for atoms and small molecules, 
the very same systems serve as diagnostic tools for the 
pulses themselves 0- The latter statement, of course, as- 
sumes that an accurate theoretical description is at hand 
for the pulsed laser-matter interaction. It is the purpose 
of this work to add further to this understanding. In par- 
ticular we are concerned with the CEPD effects, and a 
geometric interpretation of these. 

Carrier-envelope phase difference effects were stud- 
ied theoretically with emphasis on CEPD-induced spa- 
tial asymmetries in the ATI-spectrum/angular distribu- 
tion in a number of pap ers on strong- field ionization of 
atoms [1, 0, fiil U , Il2l ] and molecules [3] , strong-field 
dissociation |l4l ] , and high- harmonic generation fla j . The 



FIG. f : (Color online). The electric field E(t) /Eo, normalized 
to the peak field strength Eo and the pulse envelope f(t) as a 
function of time for two values of the CEPD, <f>. The electric 
field points in opposite directions for = and <f> = it. In ion- 
ization, the electric field will shake the electron until it gains 
enough energy to escape the Coulomb potential, and the an- 
gular distribution will depend on CEPD because the electric 
field (and the force F = —E) points in opposite directions for 
cj> — and <j> = 7T. The field parameters are r = 441 a.u. and 
lu = 0.057 a.u. (800 nm). 



asymmetries can be used to extract information about 
and ultimately to measure the CEPD. In Ref. jl6[ a spa- 
tial asymmetry in ionization with few-cycle circular po- 
larized laser pidses was observed for the first time. In 
Ref. 



17j | the generation of intense few-cycle laser pulses 



with stable CEPD was demonstrated, and a way of mea- 
suring CEPD, based on soft-X-ray radiation was pre- 
sented. In Ref. [l8| a spatial asymmetry in the ionization 
with few-cycle linear polarized laser pulses was measured, 
and the CEPD was determined with an estimated error 
of 7r/10. Asymmetries in the ionization signal combined 
with an attosecond pump pulse @ were used to measure 
directly the field of a linearly polarized few-cycle pulse 0] • 
Very recently, the detailed control over the CEPD played 
a crucial role in attosecond electron dynamics (IR Hfjl. lifj] . 
and experiments supporting the present findings have 
been reported at conferences [ilj ]. 

In this work, we present a systematic theoretical study 
of the symmetry of the response of atomic and molecular 
systems under a change of the CEPD in few-cycle pulses, 
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and we exemplify the discussion with a study of the dif- 
ferential electron momentum distribution for hydrogen 
under such pulses. 

We start out by considering an n-electron atom inter- 
acting with a few-cycle circularly polarized laser pulse 
described by the vector potential 



(1) 



x ^cos(?7 + 6 - 



^-)e x + sin(ry + . 



with f(rj) = sin (^y) the envelope, N the number of 
optical cycles, r\ = ujt — kz, uj the frequency, and k = ke z 
the wave vector. In the present case, with full inclusion of 
the spatial dependence of the field, the interaction of the 
atom with the field is obtained by the minimal coupling 
p — > p + A and the time-dependent Schrodinger equation 
reads [atomic units (a.u.) with m e — e = ao = h — 1 are 
used throughout], 
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where Hq is the field-free Hamiltonian. We are interested 
in relating this equation to an equation for the 6 = 
case. To this end, we note that the unitary operator 
D(6) = exp(jj z 0), where J z is the total angular mo- 
mentum, corresponds to a rotation of our system by an 
angle —6 around the z-axis. Since Hq is invariant un- 
der rotations around the z-axis, the transformed wave 
function = D{4>)'^ satisfies the Shrodinger equation, 

i-V =(ffo + J>(0)^,f,fi) -ffj^{4>) (3) 



nA 



^D(0)/ 2 (^)I>t(0))*', 



and using the Baker-Hausdorff lemma [22j we obtain, 
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^ = | H + ^ A(o, t, fj) ■ ft + ) 

(4) 

When we compare @ and Q , we see that a change in the 
CEPD from <f> — to 6 = 6 ' corresponds to a rotation 
of our system around the z-axis by the angle 6 ' ' ■ This 
fact is not only theoretically interesting, but also helpful 
in practical calculations since this symmetry property re- 
duces the number of computations one has to perform to 
a single one - all other results are obtained by suitable 
rotations. For instance, imagine we are interested in the 
differential ionization probability ^ ^ (<f>) for the mo- 
menta q x , q y in the polarization plane and for a general 
4> > 0. Then we calculate ^qdq = an< ^ r °t a te the 
result counterclockwise by an angle (f> to obtain d ^ P ^ q {4>) ■ 
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In the above derivation it is essential that the field-free 
Hamiltonian Hq is invariant to rotations around the z- 
axis. If this is not the case, the proof breaks down. As 
an example we look at ionization of diatomic molecules, 
or more generally linear molecules, in circularly polar- 
ized few-cycle laser pulses. If the molecule is aligned 
along the laser propagation direction, then the field-free 
Hamiltonian still has the required symmetry and the the- 
orem holds. If the molecule is not aligned along this 
axis, the system does not have the required symmetry 
and CEPD effects can not be reduced to a geometrical 
rotation. Accordingly, we may make a distinction be- 
tween rotational invariant atomic and molecular systems 
where the CEPD effects are purely geometric rotations, 
and systems which are not rotational invariant in which 
case true non- geometrical CEPD effects occur. For ex- 
ample, the results on strong-field ionization of with 
the molecule in the polarization plane [l3| belong to the 
latter category. 

As an illustration of the present findings, we consider 
ionization of atomic hydrogen, H(ls), in the strong-field 
approximation (SFA) [23|. We assume the dipole ap- 
proximation which means that rj — tot in (fTJ) and use 
that E — —d t A. The probability amplitude for direct 
ionization reads 



T fi = -i 



^ff{r,t)\E-r\Mr,t))dt, 



(5) 



where ^(r, t) is the Is ground state wave function. The 
final state f(r,t) is represented by a Coulomb- Volkov 
wave function [24j with asymptotic momentum q. The in- 
tegral in |5]| is analyzed as in [24[ and the numerical eval- 
uation is performed with Gauss-Legendre quadrature. 
By exploiting the invariance of the dot-product under ro- 
tations (E{4>) ■r*= [R z (<p)E(<p = 0)] • [R z (^)R z {-^)r\ = 
E(4> = 0) • f ', where R z (4>) is the 3x3 matrix describing 
rotations around the z axis), the SFA is readily shown to 
respond to CEPD changes like the exact theory discussed 
above. 

Once Tfi is known, a simple numerical integration over 
q z gives us the (q Xl q y ) distribution ^'^ = J\T fi \ 2 dq z . 
Figure presents the calculated distribution for various 
values of cf>. For varying cj>, the distribution rotates in 
accordance with the general theory. For 6 — 0, f Fi , % 

(IQ.r o Q u 

has a peak around (q x = 0, q y ~ 0.6) (this peak corre- 
sponds to about 14 photons above threshold), and is al- 
most symmetric around the line q x = 0. As 6 increases, 
this line is turned counterclockwise, so that has a 

peak around (q y = 0, q x ~ —0.6) for 6 = 

Some of the features can be explained by a semiclas- 
sical two-step model llj: first the electron escapes to 
the continuum at t = to with velocity v(to) = 0. Sec- 
ond, it moves like a classical particle under the influ- 
ence of the external Coulomb and laser fields. If, we 
neglect the Coulomb potential, which is justified if the 
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FIG. 4: (Color online) The {q x ,q y ) distribution for strong- 
field ionization of H(ls) by a linearly polarized field for (j> — 



0,7r/2, with I = 5.0 x 10 13 W/cm 2 



0.057 corresponding 



to a central wavelength of 800nm and three cycles, N 



The grid size is dq x = dq y 



0.02. 



FIG. 2: (Color online) The (q x ,q y ) distribution for strong- 
field ionization of H(ls) for various values of <j>, with I = 5.0 x 
10 13 W/cm 2 , to — 0.057 corresponding to a central wavelength 
of 800nm and three cycles, N = 3. The grid size is dq x = 
dq v = 0.01. 
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FIG. 3: The curves show the time-dependent vector potential 
and the lines show if* for cj> = 0, tt/2. Laser parameters are as 
in Fig. 



field is very strong, Newton's second law tells us that 
q x (t) = A x (t)+q d , q y {t) = A y (t)+q d where q d x = -A x {t ) 
and q d = — A y (to). This means that the photoelectron 
direction, after the pulse where A vanishes, is determined 
by the vector q d = (q d , q d ). For simplicity t is chosen to 

be t/2, which is the time when \E(t)\ reaches its maxi- 
mum. Figure [3] shows plots of the time-dependent vector 
potential and q d = —A(to) for various values of <j>. The 
model explains how the preferred direction(the line of 
symmetry) of the photoelectron depends on CEPD. 

Finally we consider the case of a linearly polarized laser 
field and return to the question of symmetry. The vector 
potential in the linear case, A((j),t,r) — A f(r]) sm(r] + 
+ %/2)e y , can be written as A((f>, t, r) — t) — 

A R (^,t)), where 1 L/R (<f>,t,r) = ^f( v )(cos(i 1 + + 
ir /2)e x ±sm(r] + (j)+Tr /2)e y ) describe left(+) and right (—) 
hand circularly polarized fields. This relation suggests 
that we can characterize CEPD effects for a linearly po- 



larized pulse similarly to circularly polarized pulses. Un- 
fortunately, in exact theory, it is not possible to separate 
the effects of the circularly polarized fields, since they 
do not individually commute with the field-free Hamilto- 
nian. 

Figure 4 shows two examples for a linear field correspond- 
ing to the 4> = and (j) — tt/2 cases. These results may 
be interpreted in terms of the classical model discussed 
above (q y = — A y (t ))- For <p — 0, we have the peak of 
the field at the peak of the envelope and this produces 



electrons with positive momentum, 



-Ml)' For 



<f> = 7r/2, the vector potential is antisymmetric with re- 
spect to the peak of the envelope, and the electron spec- 
trum reflects this in a forward-backward symmetry. For 
4> = tt/2, the spectrum has moved to smaller momenta 
since the two maxima of the vector potential are lower 
in this case. Notice, however, that the model can not 
explain the size of the shift. There are two reasons for 
this discrepancy. First, the model neglects the Coulomb 
potential. Second we assume that vo = 0, which is not 
necessarily correct. 

In this work, we provided an exact characterization 
of CEPD effects in the interaction between circularly 
polarized few-cycle pulses and multielectron atoms and 
molecules including full account of the spatial depen- 
dence of the field (non-dipolc effects). For systems which 
are invariant to rotations around the propagation direc- 
tion of the field, a change of the CEPD corresponds to 
a rotation of the systems. For rotationally non-invariant 
systems, true non-geometrical CEPD effects may be ob- 
served. Calculations on H(ls) quantified the effects, and 
showed that observation of the momentum distribution 
can be used to measure CEPD for a circularly polarized 
laser pulse. For linearly polarized fields, the electrons 
spectrum was interpreted in terms of a classical model. 

There is currently considerable interest in control 
of the response of atoms and small molecules to a 
strong few-cycle pulse [H, e.g., in connection with align- 
ment dependent high-harmonic generation from small 
molecules [2|| [26[. Here the harmonics provide a deli- 
cate probe of the molecular wave packet dynamics and, 



4 



in turn, steering the electronic wave packets controls the 
harmonics. The carrier-envelope-phasc difference is one 
of the essential control parameters in few-cycle pulse and 
attosecond science. 
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